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ABSTRACT 

The KATRIN experiment in Karlsruhe Germany will monitor the decay of tritium, 
which produces an electron-antineutrino. While the present upper bound for its mass is 
2 eV/c 2 , KATRIN will search down to 0.2 eV/c 2 . If the dark matter of the galaxy cluster 
Abell 1689 is modeled as degenerate isothermal fermions, the strong and weak lensing 
data may be explained by degenerate neutrinos with mass of 1.5 eV/c 2 . Strong lensing 
data beyond 275 kpc put tension on the standard cold dark matter interpretation. In 
the most natural scenario, the electron antineutrino will have a mass of 1.5 eV/c 2 , a 
value that will be tested in KATRIN. 

Subject headings: Introduction, NFW profile for cold dark matter, Isothermal neutrino 
mass models, Applications of the NFW and isothermal neutrino models, Fit to the most 
recent data sets, Conclusion 

1. Introduction 

The neutrino sector of the Standard Model of elementary particles is one of today's most intense 
research fields (Kusenko, 2010; Altarelli and Feruglio 2010; Avignone et al. 2008). While neutrino 
oscillations give information of differences in mass-squared values, no absolute neutrino mass is 
known. The masses may be of order 0.05 eV/c 2 as set by the oscillations, or much larger. While 
a lot is known about standard (active) neutrinos (left-handed neutrinos, right-handed antineutri- 
nos), hardly anything is known about their sterile partners (right-handed neutrinos, left-handed 
antineutrinos) . It is even not established that the masses of active and sterile components are close 
to each other; in the popular seesaw mechanism, they differ by orders of magnitude (Lesgourges 
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and Pastor 2006). Still, next to WIMPs (Bertone 2010), sterile neutrinos may constitute the dark 
matter particles (Feng 2010). It was recently realized that lensing data allow that active plus sterile 
neutrinos of 1.5 eV/c 2 provide the dark matter (Nieuwenhuizen et al. 2009). Neutrinos are too 
light to trigger large scale structure formation (galaxy formation), and are therefore often ruled 
out as dark matter candidates. However, this may end up as a misconception, since structure 
formation can actually already be explained by the gravitational hydrodynamics of baryons alone 
(Nieuwenhuizen et al. 2009). 

The most effective way to determine the mass of an antineutrino is tritium decay, when one of 
the two neutrons in its nucleus decays into a proton, an electron and an electron-antineutrino. The 
electron energy can be detected. Analyzing the number of events (that is, the Kurie function) near 
the maximal energy provides information about the antineutrino mass. Depending on statistical 
assumptions in the analysis of the non-detection in the Mainz and Troitsk experiments, the present 
upper bound is 2.2 eV/c 2 (Lobashev et al. 1999), 2.3 eV/c 2 (Kraus et al. 2005) or 2.0 eV/c 2 
(Amsler et al. 2008). 

The KArlsruhe TRItium Neutrino experiments (KATRIN) involve an international consortium 
of 16 laboratories in Europe and the USA, and are intended to search the mass down to 0.2 eV/c 2 
(Often and Weinheimer 2009). Since the square mass is tested, this means to go down to 0.04 
eV 2 /c 4 , an improvement by two orders of magnitude with respect to Mainz-Troitsk. 

Galaxy clusters are known to produce strong lensing effects. A galaxy lying behind it can 
be seen as an arc, ideally as a complete Einstein ring. This has been observed, see Fig. 1 which 
presents a double Einstein ring. The very similar phenomenon may be seen already in lensing by a 
glass ball, see figure 2. This so-called strong lensing is effective near the center of the cluster, say 
up to a radius of 300 kpc. Further out, there is still the so-called weak lensing. Indeed, since every 
galaxy will be slightly distorted, a statistical effect occurs for a "local group of galaxies ("local on 
the 2D sky, as observed by us), which then also provides information on the mass distribution of 
the lensing cluster. Models for structure of the cluster can be tested against these strong and weak 
lensing data. The clusters also have an amount of 10-15% of matter in X-ray gas, typically much 
more than they have galaxy mass. 

The galaxy cluster Abell 1689 (A1689) is one of the most studied clusters. It is rather heavy, 
a few times 10 15 solar masses, and has a large lensing arc at 150 h^Q kpc, which is observed at 45" 
from the center, the largest Einstein radius observed to date. (The Hubble constant is Hq = Zi7o70 
km/s Mpc, with hjo close to 1.0). There are good data for strong lensing (SL) and weak lensing 
(WL) and the X-ray gas, so it is an ideal testing ground for models. The cluster has a small intruder 
clump in the North East sector, which can be masked out if needed. Assuming spherical symmetry, 
we shall consider two models for the matter distribution in the cluster. 
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Fig. 1. — Strong lensing by a foreground galaxy photographed by the Hubble Space Telescope HST. 
Two background galaxies are mapped as two almost complete Einstein rings. In the enlargement 
on the right, the foreground galaxy is masked out. 




Fig. 2. — Strong lensing by a glass ball at home. The central square of the background paper is 
distorted. Parallel lines are mapped as semi-arcs, pairs of which make up the ellipses. 

2. NFW profile for cold dark matter 

Cold dark matter (CDM) is mostly thought to be related to particles heavier than the pro- 
ton, like the lightest supersymmetric (SUSY) particle called neutralino. However, the assumption 
of CDM is subject to criticism. The so-called R-symmetry assumption implies a very long life- 
time. The electron and proton have spin, are both charged and light, so they "have nowhere to 
go"; but if the R-symmetry is broken, the neutralino would quickly decay to electrons, protons, 
neutrinos and photons. In practice, the purported CDM particle has never been found, despite 
decades of searching. An incomplete list of acronyms of the over 60 search programs is, not men- 
tioning updates: ADMX, AMANDA, ANAIS, ANTARES, ArDM, ATIC, BPRS, CANGAROO, 



- 4 - 



CAST, CDMS, CLEAN, CMSSM, CoGeNT, COUPP, CRESST, CUORE, CYGNUS, DAMA/Nal, 
DAMA/LIBRA, DAMIC, DEAP, DMTPC, DRIFT, EDELWEISS, EGRET, ELEGANTS, FERMI- 
LAT, GENIUS, GERDA, GEDEON, HDMS, ICECUBE, IGEX, KIMS, LEP, LHC, LIBRA, LUX, 
MAGIC, MAJORANA, MIMAC, NAIAD, NEWAGE, ORPHEUS, PAMELA, PICASSO, ROSE- 
BUD, SIGN, SIMPLE, SuperCDMS, SUPER-K, TEXONO, UKDM, VERITAS, WArP, Whipple, 
XENON, XMASS, ZEPLIN. Upcoming are: DARWIN, DEAP-3600, EURECA, GEODM, Mini- 
CLEAN, SuperCDMS-SNOLab. (Notice that J, Q and Y are not covered.) Also in recent years 
the searching remains without detection, in particular the XENON 100 search rules out previous 
detection claims at DAMA, CDMS and Cogent (Aprile et al. 2010), a rinding confirmed by the 
CDMS collaboration (Akerbib et al. 2010). While it is already acknowledged that time for su- 
persymmetry is running out (Bertone 2010), the CMS collaboration at the Large Hadron Collider 
LHC reports, very recently, that within error limits it fails to detect more SUSY-like events than 
expected from the non-SUSY background (CMS 2011). So, though the SUSY assumption and 
the possibly resulting neutralino CDM are certainly not ruled out, these long trusted options are 
definitely on the defensive at the moment. 

From the practical side in cosmology, the CDM assumption works often so well that it is 
presently viewed as the standard approach. This applies in particular for the cosmic microwave 
background. As to cosmic structures, numerical simulations have shown that the CDM density 
can be described fairly well by an NFW profile (Navarro, Frenk and White, 1997), which has the 
benefit of having a simple form. It is parametrized by a characteristic mass density po an d a length 
scale r s , 



The average density, that is, the integrated mass up to R divided by the volume, then equals 



The cosmic critical density is p c = 3H$/8ttG = 9.25 • 10- 27 kg/m 3 according to the WMAP7 
value for the Hubble constant, H = 70.2 ± 1.4 km/s/Mpc (Komatsu et al. 2010). At the redshift 
z = 0.183 of A1689 the Hubble constant reads H z = E z H with 



Pnfw(0 = 



Por 3 s 



(1) 



r(r + r s ) 2 ' 




(2) 



E z = [n M (l + zf + (1 - n M - O0(l + z? + «a] 1/2 = 1-078, 



(3) 



(we adopt the flat cosmology U\ = 0.75, Qm = 0.25), so the critical density becomes p c (z) = p c E z . 
The radius Rjy where /?nfw(-^/v) equals N p c (z) is said to be the radius where the over-density 
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equals N. It defines the concentration parameter cat = R^/r s . Customary values are N = 200 
corresponding to C200 and "virial" overdensity N = 117 at the redshift of A1689, with its c v ; r = c\yj. 
The correspondence involved here is 

P° = 3^v , 4) 

N Pc (z) log(l + c N ) - c N /(l + c N ) ' 

with the monotonically increasing right hand side yielding a unique cjy. 



3. Isothermal neutrino mass models 



It is long believed that light neutrinos cannot make up the dark matter, because they cannot 
assist in structure formation of galaxies and clusters of galaxies. Indeed, the Jeans mechanism in hy- 
drodynamics can only explain globular star clusters of, say, 600,000 solar masses. But gravitational 
hydrodynamics offers new mechanisms for cosmological structure formation, with an important role 
for viscosity. We shall not dwell into it here, as it is discussed at length in other papers in this 
volume by Gibson and Schild (2011) and Gibson (2011). This approach also answers the paradox 
why the 3 He abundance in the Galaxy is often observed at primordial level, independent of the 
the distance to the center and thus independent of the local amount of metallicity (Nieuwenhuizen, 
2011). For the present work it suffices to notice that neutrinos need not be ruled out by arguments 
of structure formation and hence qualify "again" for dark matter candidates. This holds the more 
since the early runs of the Large Hadron Collider at the CERN did not detect supersymmetry 
(CMS 2011). 

One may thus consider the cluster dark matter as due to degenerate fermions. The violent 
relaxation mechanism has explained the success of isothermal models in gravitation (Lynden Bell 
1967), so it is natural to consider isothermal dark matter fermions. 

The gravitational potential caused by mass elements at r' / r is ip(r) = —G f d 3 r'p(r') /|r — r'| 
+ G J d 3 r / p(r / )/r / , where we choose the zero level such that <p(0) = 0. It satisfies the Poisson 
equation V 2 (^(r) = 47rGp(r), which in case of spherical symmetry reads 



(p"(r) + -tp'(r) = 47rG/9(r), (5) 

In the galaxy cluster there will be three relevant contributions to the matter density, p = Pu+PG+Pg-, 
with partial mass densities of dark matter ('V"), Galaxies (G) and X-ray gas (g), respectively. We 
model the dark matter as isothermal fermions with mass m, degeneracy g and chemical potential 
p = aksT. The mass density then involves a momentum integral over a Fermi-Dirac distribution, 



Pv{ r ) = gm 



I 



d 3 p 1 
(2irh) 3 exp{[p 2 /2m + mtp(r) - p]/k B T} + 1 ' 



(6) 
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The thermal length is At = fry / 27r / 'rnksT '. 

Since the mass in galaxies is a few percent of the total, mainly coming from the dominant 
central one, we assume for the mass density from galaxies simply the isothermal distribution pc(r) = 
gmX^ 3 exp[«G — /3cmip(r)] (Nieuwenhuizen 2009). In the galaxy cluster Abell 1689 the mass density 
of the gas p g follows from the electron density n e as p g = 1.11 m^n e with mjv the nucleon mass. 
The data for n e densely obtained by us from X-ray observations in the interval 10 kpc < r < 1 
Mpc, are well approximated by the exponent of a parabola in log r (Nieuwenhuizen and Morandi, 
2011) 

n e (r) = exp( -3.63 + 0.915 log - 0.218 log 2 \ hjjV" 3 . (7) 

V Kpc Kpc / 

Beyond 1 Mpc one may either continue this shape, or replace it for r > r- lso = 807 kpc by the 
isothermal tail (Nieuwenhuizen and Morandi, 2011) 

""=2^' ^ = 578.5^. (8) 

From strong lensing data one determines the 2D mass density (integrated mass along the line 
of sight) S(r) = dzp (^Vr 2 + z 2 ^ . The Poisson equation allows to express this as 



, > 1 f°° , cosh 2s 



(p [r cosh s) 



cosh 2 s 



(9) 



which is numerically more stable since (p' varies less than p. In weak lensing one determines the 
shear as 



where the "critical projected mass density" S c is a constant. The mass inside a cylinder of radius r 
around the cluster center, M2_o(r) = 2ir dr' r'S(r'), defines the average of S over a disk of radius 
r, S(r) = M2rj(r)/7rr 2 . The Poisson equation allows to express this as (Nieuwenhuizen 2009) 

— 1 f°° 

£(r) = — / dst^rcoshs). (11) 
k(j Jo 



4. Applications of the NFW and isothermal neutrino models 



Cowsik and McClelland (1973) are the first to model the dark matter of the Coma cluster 
as an isothermal sphere of neutrinos. As for neutron stars and white dwarfs, this determines the 
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neutrino mass, and yields about 2 eV/c 2 . Treumann et al. (2000) considered 2 eV/c 2 neutrinos 
next to CDM and X-ray gas, all isothermal, and derive the density profiles for clusters such as 
Coma. Nieuwenhuizen (2009) applies to clusters an isothermal fermion model for a single type of 
dark matter, cold or not, and also isothermal for the galaxies and the X-ray gas. Fitting to weak 
lensing data of the Abell 1689 cluster works well and yields as the mass (I2/5) 1 / 4 1.455 ± 0.030 /4o 2 
cV/c 2 , where g is the fermion degeneracy. The only known candidate for such a low mass is the 
neutrino, for which g = 3*2*2 = 12 modes are available, due to 3 families of (anti)particles 
with spin ^, because helicity is not conserved in the gravitational field of the cluster. The mass 
lies below the experimental upper bound of 2.0 eV/c 2 but within the KATRIN window. With this 
mass active neutrinos constitute about 10% dark matter. Some 20% neutrino dark matter may 
arise provided the sterile partners have been created, too, in the early Universe. This may occur in 
the presence of a me V- valued Majorana mass matrix, in the time window between the freeze out 
of quarks and the freeze out of electrons and positrons (Nieuwenhuizen 2009). 

Umetsu and Broadhurst (2008) combine weak lensing data with reddening of background 
galaxies out to 3 Mpc, presenting the best (extended) weak lensing data to date. Combining 
them with certain strong lensing data they fit an NFW profile and find a concentration parameter 
c 200 = 13.5, well exceeding the expected value around 6-7. Nieuwenhuizen and Morandi (2011) 
employ the Umetsu and Broadhurst (2008) weak lensing data but use better strong lensing data 
of Limousin et al. (2007) and arrive at the more expected C200 = 6.8. However, they point at a 
systematic deviation in the weak lensing fit, related to the fact that weak lensing alone would give 
a value of around 7.6 (Limousin et al. 2007). They also point at a mild tendency of the strong 
lensing data to fall well below NFW beyond 300 kpc. 

Coe et al. (2010) present improved strong lensing data by employing a non-parametirc ap- 
proach. Combining these strong lensing data up to 220 kpc with the weak lensing data of Umetsu 
and Broadhurst (2008) beyond 220 kpc, they fit to an NFW profile and get r s = 338 kpc, C200 = 7.6 
(or c vir = 9.6), leading to r M 200 = IA^IQ^K^Mq or M vir = 2.Ot^lO 15 /if o 1 M . 

Nieuwenhuizen and Morandi (2011) also apply the isothermal neutrino model to their A1689 
data. The approach again works well, with slightly modified parameters and neutrino mass 1.499 ± 
0.039 h^^ 2 eV/c 2 . This theory supports a fast decay of strong lensing beyond 300 kpc. The reason 
is that its neutrino dark matter would be localized up to, say, 300 kpc (Nieuwenhuizen 2009), with 
the X-ray gas more concentrated in the outskirts. As seen in Fig 4 for a slightly different data set, 
in this model strong lensing decays fast beyond the region where DM is localized. Weak lensing, 
on the other hand, does not exhibit this, as it senses a slow decay, namely an isothermal 1/r decay 
plus an 1/r 2 correction (Nieuwenhuizen and Morandi, 2011), see also Fig. 5 for weak lensing in the 
data set considered below. 

The case of neutrino dark matter in clusters is also considered by Sanders (2007). While he 
assumes modifications of Newton gravity, we shall stay within the Newtonian scope. 
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Fig. 3. — Reproduction of Figure 2 of Coe et al. (2010) by permission of the AAS. Blue: Strong 
lensing data of the cluster A1689 obtained with their analysis. Black: binned data. Green: data 
by Broahurst et al. (2005). Red: The NFW model fits well, but it does not support the fast fall 
off beyond 275 kpc. 

5. Fit to the most recent data sets 

We now consider the best data sets to date, the strong lensing of Coe et al. (2010), combined 
with the weak lensing data of Umetsu and Broadhurst (2008) and the X-ray data of Nieuwenhuizen 
and Morandi (2011). Whereas Coe et al. (2010) used the strong lensing data up to 220 kpc and 
the weak lensing data beyond this, we shall now follow our 2011 work and use all strong lensing 
and weak lensing data, taking the gas data as input. We wish to minimize x 2 = Xsl + Xwl f° r * ne 
20 strong lensing data points and 11 weak lensing points. 

Let us first discuss the NFW case, which involves 2 fit parameters, leaving a minimization over 
20+11-2=29 non-trivial inputs. (Not having the correlation matrix between the data points of SL 
or WL, we perform x 2 fits, here and further on.) The fit leads to Xsl = ^1.4 and Xwl = 14.6. 
The profile parameters are r s = 291 kpc and C200 = 8.8, or the virial value c v i r = 11.0. With its 
X 2 /29 = 1.24, this fit is satisfactory. But, as seen in Fig 4, the strong lensing strongly falls down 
beyond r s , while the NFW profile continues gradually. So if the fall off is statistically relevant, it 
will pose a problem for the NFW model. 

One can also make a fit to the isothermal fermion model assuming an isothermal distribution 
for the density of matter in galaxies, pair) = gmX^ 3 exp[ac — (3G m( p( r )] (Nieuwenhuizen 2009). 
The fit involves five parameters, the best values being 
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a = 29.26 ±0.43, ?J = 4092 ± 275^, aj, = 564 ±9—, 
Aj, ur s 

« G = 7.44 ±0.11, fc-!™^, (12, 

where a v = ^/ksT/m is a DM velocity dispersion for which At = h\/2ir/m<j u . The result for fie 
differs the most with respect to fits (Nieuwenhuizen and Morandi, 20011), but it still is close to 
1/mal, which is the case of virial equilibrium of galaxies and DM (Nieuwenhuizen 2009). The mass 
comes out as 



/12\ 1/4 

m=[ — \ 1.55eV±0.04eV/c 2 . (13) 

From now on, we take g = 12; in general the masses should be multiplied by the factor (12/g) 1 / 4 . 

The strong lensing fit, shown in Fig. 4, corresponds to Xsl = 57.9. Its match below 100 kpc 
could be improved by a different modeling of the galaxy density, which would also allow to improve 
the fit between 100 and 300 kpc, so this could easily bring down the present large value of x!l- 
The weak lensing fit is shown in Fig. 5. The outlying 8'th weak lensing point alone contributes an 
amount 10.0 to Xwl = 18.1, the 10 other points bringing 8.1, so that weak lensing is generally fit 
quite well. In total, the present fit thus gives the value \ 2 1 (31 — 5) = 2.92, larger that for NFW. We 
shall not attempt to improve the fit, since this will have little impact on the value of the neutrino 
mass. Moreover, it is seen in Fig. 4 that the neutrino model predicts a fast decay of strong lensing 
beyond 275 kpc, while NFW does not. So if this decay of the Coe et al. (2010) data in Fig. 3 is 
statistically relevant, it supports the neutrino model and disfavors NFW. In this regard one may 
determine till which radius NFW predicts lensing arcs, and compare this with the observations. 



6. Conclusion 



The lensing data of the galaxy cluster Abell 1689 may be explained by isothermal fermions. The 
best case is neutrinos of nearly the same mass. Combining with previous (but less good) data sets 
on the same cluster we infer a mass is m u = 1.50 ±0.05 (sys) ±0.03 (stat) eV/c 2 . The three species 
of active neutrinos have been created in thermal equilibrium and they are predicted have total mass 
of M v = 2>m v ps 4.5 eV/c 2 . With the general relation Q v = M u /(46hj ) the mass density in active 
neutrinos would be f^ ct ps 10% of the critical mass density of the Universe. This is much more than 
the 1.2-3.0% allowed by the standard ACDM approach under various conditions and combinations 
with other data sets (Komatsu et al. 2010). In our philosophy neutrinos should explain most, if 
not all, cluster dark matter, Q, u = 20 — 25%. Under reasonable conditions, another 10% can have 
arisen from the creation of three thermal sterile neutrinos species of 1.5 eV/c 2 (Nieuwenhuizen, 
2009). If there are more sterile than active partners that would explain even more dark matter. 
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Fig. 4. — Dots: log-log plot of the binned strong lensing data of Coe et al. (2010), see Fig. 3. 
Dashed line: The NFW model fits well, but fails to explain the fast fall off beyond 275 kpc. Full line: 
Isothermal dark fermions with the measured gas and isothermal modeling of galaxies fit reasonably. 
Dotted lines: the separate contributions from (neutrino) dark matter, Galaxies and X-ray gas. 
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Fig. 5. — Weak lensing data of Umetsu and Broadhurst (2008). Full line: isothermal neutrinos. 
Dotted line: its asymptotic shape. Dashed line: NFW profile of cold dark matter. 

The 1.5 eV/c 2 mass lies in the range covered by the KATRIN experiment, searching the mass 
of the electron antineutrino, that is expected to start taking data in 2011 (Weinheimer 2009). A 
discovery around 1.5 eV/c 2 would clearly prove our thesis. But non-detection would not necessarily 
rule out the case of nearly all non-baryonic dark matter in neutrinos. What KATRIN measures is, 
in leading approximation, the combination m 2 e = ^ \ U e i\ 2 mf, where U e i is the amplitude of the 
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neutrino mass eigenstate \i>i) in the electron antineutrino wavefunction, \u e ) = Yli^eiWi), with i 
running from 1 to 3 + N s with N s the number of sterile neutrinos, and rrij the mass of the eigenstate 
\vi) (Otten and Weinheimer, 2009). The LSDN (Aguilar et al. 2001) and MiniBoone data for short 
baseline neutrino oscillation experiments (Aguilar et al. 2010) can be explained by assuming one 
or more sterile neutrinos with mass of order 1 eV/c 2 , with the solar and atmospheric neutrinos 
involved with much lighter masses. Therefore a conspiracy may take place, the heavier mass 
eigenstates having small matrix elements U e i, making m 2 e too small to be detected at KATRIN, 
while dark matter may still come from neutrinos, namely the sterile ones, with masses of order of 
eV/c 2 . A more interesting scenario is that the masses mi are well enough separated for KATRIN to 
observe their effect one-by-one. (Otten and Weinheimer, 2009). Such questions and many related 
ones are beyond the scope of the present paper, but see Nieuwenhuizen (2011). Clearly, the (sterile) 
neutrino sector of the standard model of elementary particles is an exciting field of research for the 
near future. 

As generally expected, the lensing data in A1689 can also be explained by an NFW profile. 
One may wonder whether there is a discriminating feature between NFW and isothermal neutrinos, 
that is, between CDM and i/DM. There indeed is a candidate in the fast fall off of strong lensing 
beyond 270 kpc, first pointed at by Nieuwenhuizen and Morandi (2011) and more convincingly 
present in the data of Coe et al. (2010), see Figs. 3, 4. NFW would predict a smooth density 
profile and smooth continuation of strong lensing data. The observations conflict with this and give 
better support for the isothermal neutrino model, where this fall off is related to the fast fall off of 
neutrino density beyond, say 300 kpc. The fact that neutrinos are more localized than the X-ray 
gas also explains the "missing baryon paradox" for this cluster, because the X-ray gas is generally 
more concentrated in the outskirts, as is concluded from observations on a set of clusters (Rasheed 
et al 2010). One may investigate whether an NFW profile predicts beyond 300 kpc more lensing 
arcs than the few observed ones. The answer to this question may touch the validity of the NFW 
model and its assumed cold dark matter. 

A.M. acknowledges support by Israel Science Foundation grant 823/09. 
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